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HYCOM  2.2  (I) 


Maintain  all  features  of  HYCOM  2.1 

o  Orthogonal  curvilinear  grids 

o  Can  emulate  Z  or  Sigma  or  Sigma-Z  models 

o  It  is  “Arbitrary  Lagrangian-Eulerian”,  see: 
Adcroft  and  Hallberg,  O.  Modelling  1 1  224-233. 

o  Explicit  support  for  1-D  and  2-D  domains 

o  KPP  or  Kraus-Turner  or 
Mellor-Yamada  2.5  or  Price-Weller-Pinkel 

o  Rivers  as  bogused  surface  precipitation 

o  Multiple  tracers 

o  Off-line  one-way  nesting 

o  Scalability  via  OpenMP  or  MPI  or  both 

o  Bit-for-bit  multi-cpu  reproducibility 


HYCOM  2.2  (Ha) 

Alternative  scalar  advection  techniques 

o  Provided  by  Mohamed  Iskandarani 
o  Donor  Cell,  FCT  (2nd  and  4th  order),  MPDATA 
o  FCT2  replaces  MPDATA  as  standard  scheme 

Vertical  coordinate  changes 

o  Vertical  remapping  uses  PLM  for  fixed 
coordinate  layers 

o  Thin  deep  iso-pycnal  layers 

o  Stability  from  locally  referenced  potential  density 

o  Spatially  varying  layer  target  densities 

o  Different  isopycnal  layers  in  semi-enclosed  seas 


HYCOM  2.2  (lib) 

Equation  of  state  that  is  quadratic  in  salinity 

o  HYCOM  must  “invert”  the  equation  of  state 

o  tofsig(r,s)  and  sofsig(r,t) 

o  Traditional  version  is  cubic  in  T  and  linear  in  S 

o  Finding  the  root  of  a  cubic  is  expensive, 
but  exact 

o  Linear  in  S  is  not  accuarate  at  low  salinity 
o  Optional  version  is  cubic  in  T  and  quadratic  in  S 
o  Coefficents  provided  by  Shan  Sun 
o  More  accurate  at  low  salinity 
•  Rivers,  Black  Sea,  Caspian  Sea 
o  Not  much  more  expensive 


HYCOM  2.2  (lie) 

Special  halo  exchange  for  tripole  global  grid 

o  Arctic  dipole  patch  on  standard  Mercator  globe 
o  Logically  rectangular  domain 
o  Two  halves  of  top  edge  “fold”  together 
o  V-velocity  changes  sign  across  the  fold 

Improved  thermobaricity 

o  No  single  reference  state  is  approriate 
for  the  global  ocean 

o  Hallberg,  Ocean  Modelling,  8,  279-300 

o  Use  a  linear  combination  of  pressure  gradients 
from  two  out  of  three  reference  states 

o  Atlantic  (3°C,  35.0  psu) 

o  Arctic/Antarctic  (0°C,  34.5  psu) 
o  Mediterranean  (13°C,  38.5  psu) 
o  Most  locations  use  just  one  reference  state 

o  Linear  combinations  allow  smooth  transition 
between  states 

•  Do  this  in  shallow  water  if  possible 

o  In  deep  water,  constrain  the  T&S  used  for 
thermobaricity  to  be  close  to  the  reference  state 


1/12°  GLOBAL  THERMOBARIC  REFERENCE 

STATE 


TBARIC  MAP  (1=Arctic,2=Atlantic,3=Med.) 


HYCOM  2.2  (Ilia) 

Mixed  layer  changes 

o  GISS  mixed  layer  model 
o  Provided  by  Armando  Howard 
o  KPP  bottom  boundary  layer 
o  Provided  by  George  Halliwell 
o  KPP  tuning 

Atmospheric  forcing  changes 

o  Option  to  input  ustar  fields 
o  Best  option  for  monthly  forcing 
o  Otherwise  calculated  from  wind  stress  or  speed 
o  Can  relax  to  observed  SST  fields 
o  Improved  COARE  3.0  bulk  exchange  coefficients 
o  Black-body  correction  to  longwave  flux 
o  Climatological  heat  flux  offset,  Qc 

Q  =  ( Qsw  —  Qiw)  +  ( Qi  +  Qs)  +  Qc 
o  Qc  is  constant  in  time 

•  Typically  based  on  the  model’s 
climatological  SST  error,  times 
(say)  -45  Wm~2/°C 


deg  C 


Qc  ADDED  AFTER  FIVE  YEARS 
GLOBAL  MEAN  SST  and  SSS 


GLBa0.72,  Basin-Wide  Mean  SST 


GLBa0.72,  Basin-Wide  Mean  SSS 


LULU 


Qc  ADDED  AFTER  FIVE  YEARS 
GLOBAL  MEAN  SST  and  SSH 


GLBa0.72,  Basin-Wide  Mean  SST 


GLBa0.72,  Basin-Wide  Mean  SSH 


HYCOM  2.2  (Illb) 

Improved  support  for  rivers 

o  Still  bogused  surface  precipitation 
o  High  frequency  inter-annual  river  flow  allowed 
o  Add  it  to  atmospheric  precip,  off-line 
o  Instead  of  monthly  climatology, 
or  in-addition  to  it  (flow  anomalies) 

o  Better  control  of  low  salinity  profiles 
o  Option  for  mass  (vs  salinity)  flux 
o  Equation  of  state  that  is  quadratic  in  salinity 

Tidal  forcing 
o  Provided  by  NCEP 

o  Body  forcing  and  open  boundary  forcing 
o  Boundry  forcing  currently  for  “Flather”  ports 
o  Extend  it  to  Browning-Kreiss  ports  and  nesting 


HYCOM  2.2  (IIIc) 

New  diagnostics  within  HYCOM 

o  Time-averaged  fields  (in  archive  files) 
o  Synthetic  instrumentation 
o  Provided  by  George  Halliwell 
o  3-D  particle  tracking 
o  surface  and  constant  depth  drifters 
o  isopycnic  drifters 
o  fixed  instruments  and  moorings 


HYCOM  2.2  (Illd) 

Finer  control  over  energy  loan  ice  model 

o  Melting  point  can  be  linear  in  salinity 
o  Set  ice  minimum  and  maximum  thickness 
o  Set  ice  vertical  temperature  gradient 
o  Or  get  ice  surface  temperature  from  Ta 
o  Made  compatible  with  coupled  sea-ice  approach 

Two-way  coupling  to  LANUs  CICE  sea  ice  model 

o  HYCOM  exports: 
o  SST,  SSS,  SSH 
o  Surface  Currents 
o  Available  Freeze/Melt  Heat  Flux 
o  CICE  exports: 
o  Ice  Concentration 
o  Ice-Ocean  Stress 

o  Actual  Freeze/Melt  Heat/Salt/Mass  Flux 
o  Solar  Radiation  at  Ice  Base 


HYCOM  AND  ESMF 

Earth  System  Modeling  Framework 
http  ://www.esmf.ucar.edu/ 

o  Superstructure  couples  components 
o  Air/Ocean/Ice/Land 
o  Asynchronous  I/O  component 

•  Run  “concurent”  with  model  components 

o  Infrastructure  provides  data  structures 
and  utilities  for  building  scalable  models 

Added  a  superstructure  “cap”  to  HYCOM 

o  Simplifies  coupled  systems 
o  HYCOM  coupled  to  LANL  CICE  sea-ice 
o  Convert  atmospheric  field  processing  and  the 
energy-loan  ice  model  into  ESMF  components 

o  Use  ESMF  for  (user-level  asynchronous)  I/O 

o  Interoperate  with  other  ESMF  compliant  ocean 
models  (e.g.  Poseidon,  MITgcm,  MOM4) 

This  initial  ESMF  support  is  optional 

o  TYPE=esmf  instead  of  TYPE=mpi 
ESMF  will  eventially  required  to  run  HYCOM 

o  HYCOM  version  3.0 


HYCOM  2.2  (IV) 

Climatological  nesting  now  allowed 

o  Start  from  monthly  mean  outer  model  archive  files 
o  Allows  nested  runs  longer  than  the  outer  run 
o  But  with  less  accurate  boundary  state 
o  Probably  only  suitable  for  regional  nests 

Nesting  no  longer  requires  co-located  grids 

o  General  archive  to  archive  horizontal 
interpolation  (curvilinear) 

Hybrid  to  fixed  vertical  grid  remapper 

o  Allows  fixed-coordinate  nests  inside  hybrid 
coordinate  outer  domains 

o  HYCOM  to  (fixed-grid)  HYCOM 
o  HYCOM  to  NCOM 


HYCOM  2.2  (V) 

Enhanced  hycomproc  and  fieldproc 

o  NCAR-graphics  based 
o  Many  more  color  palette  options 
o  Can  read  in  an  arbitrary  palette 
o  Mark  locations,  and  draw  tracks,  on  plot 
o  Plot  diffusion  coefficients  and  tracers  (hycomproc) 
o  Overlay  vector  and  line-contours  (fieldproc) 

Added  fieldcell 

o  Like  fieldproc,  but  for  cell-array  (vs  contouring) 
o  Mark  locations  and  draw  tracks 
o  Overlay  line-contours 
o  Uses  NCAR’s  map  projections 

o  Typically  much  faster  than  fieldproc, 
but  can  leave  unfilled  cells 

o  Option  to  increase  resolution  of  input  (bi-linear 
interpolation) 


ERA-15  SST:  Feb 


SSS  (hour  0) 


HYCOM  2.2  (VI) 

Diagnostic  fields  to  netCDF  and  other  file  formats 

o  Archive  fields  in  layer  space 
o  On  p-grid  (interpolated  velocity) 
o  3-D  archive  fields  interpolated  to  z-space 
o  On  p-grid,  or 

o  Sampled  along  arbitrary  tracks 
o  3-D  archive  fields  sampled  on  iso-therms 
o  Meridional  stream-function  from  (mean)  3-D  archive 
o  In  logical  array  space  (rectilinear  grids) 
o  Binned  to  latitude  bands  (curvilinear  grids) 
o  Atmospheric  forcing  input  fields 

o  Time  axis  depends  on  “.b”  file  format 
o  Any  “.a”  file  with  the  right  “.b”  file  structure 
can  be  converted  to  netCDF 

o  Fields  binned  into  lon-lat  cells 


HYCOM  CURVILINEAR  GRIDS  and  NetCDF 
Most  basin-scale  cases  use  a  Mercator  grid 

o  1-D  lat  &  Ion  axes  (rectilinear) 
o  Handled  well  by  many  netCDF  packages 

Global  HYCOM’s  Arctic  patch  grid  is  curvilinear 

HYCOM  netCDF  use  the  CF-1 .0  conventions,  which 
support  curvilinear  grids 

o  If  latitude  and  longitude  are  2-D  grids 
o  1  -D  axes  are  array  indexes 
o  Longitude  and  latitude  arrays  are  also  in  the 
file  and  identified  as  alternative  coordinates 

Most  netCDF  packages  are  not  CF-1 .0  aware 

o  Bin  into  uniform  lon-lat  cells  off-line 

o  Interpolate  to  a  1-D  latitude  and  longitude 
grid  off-line 

o  General  archive  to  archive  horizontal 
interpolation 

Archive  to  archive  remapper  can  also  be  used  for 
standard  (non-native)  grids 

o  MERSEA  grid  is  uniform  1/8° 
o  AOMIP  grid  is  rotated  uniform  1/2° 


CANDIDATE  FEATURES  FOR  HYCOM  2.3 


Stable-code  vs  new  features 

o  Released  code-base  has  to  be  tested  and  stable 
o  New  features  can  be  a  significant  improvement 
o  Will  add  interim  releases  to  web  page 
o  Features  may  be  removed  in  next  released  code 

Fully  region-independent 

o  Compile  once,  run  on  any  region  and  any 
number  of  processors 

o  Run-time  memory  allocation 
o  Might  reduce  performance  (fewer  compiler 
optimizations  available) 

o  Needed  for  full  ESMF  compliance 
Wetting  and  Drying 
Diurnal  heat  flux  cycle 
Wind  drag  coefficient  based  on  model  SST 
Support  for  HYCOM  in  CCSM 
Enhanced  support  for  ESMF 


